Search for: All records

Abstract Near the ends of their lives, supernova remnants (SNRs) enter a “radiative phase,” when efficient cooling of the postshock gas slows expansion. Understanding SNR evolution at this stage is crucial for estimating feedback in galaxies, as SNRs are expected to release energy and momentum into the interstellar medium near the ends of their lives. A standard prediction of SNR evolutionary models is that the onset of the radiative stage precipitates the formation of a dense shell behind the forward shock. In Paper I, we showed that such shell formation yields detectable nonthermal radiation from radio toγ-rays, most notably emission brightening by nearly 2 orders of magnitude. However, there remains no observational evidence for such brightening, suggesting that this standard prediction needs to be investigated. In this paper, we perform magnetohydrodynamic simulations of SNR evolution through the radiative stage, including cosmic rays (CRs) and magnetic fields to assess their dynamical roles. We find that both sources of nonthermal pressure impede shell formation, reducing shell densities by a factor of a few to more than an order of magnitude. We also use a self-consistent model of particle acceleration to estimate the nonthermal emission from these modified SNRs and demonstrate that, for reasonable CR acceleration efficiencies and magnetic field strengths, the nonthermal signatures of shell formation can all but disappear. We therefore conclude that the absence of observational signatures of shell formation represents strong evidence that nonthermal pressures from CRs and magnetic fields play a critical dynamical role in late-stage SNR evolution. 

Abstract We present the results from our extensive hard-to-soft X-ray (NuSTAR, Swift-XRT, XMM-Newton, Chandra) and meter-to-millimeter-wave radio (Giant Metrewave Radio Telescope, Very Large Array, NOEMA) monitoring campaign of the very nearby (d = 6.9 Mpc) Type II supernova (SN) 2023ixf spanning ≈4–165 days post-explosion. This unprecedented data set enables inferences on the explosion’s circumstellar medium (CSM) density and geometry. In particular, we find that the luminous X-ray emission is well modeled by thermal free–free radiation from the forward shock with rapidly decreasing photoelectric absorption with time. The radio spectrum is dominated by synchrotron radiation from the same shock. Similar to the X-rays, the level of free–free absorption affecting the radio spectrum rapidly decreases with time as a consequence of the shock propagation into the dense CSM. While the X-ray and the radio modeling independently support the presence of a dense medium corresponding to an effective mass-loss rate $\dot{M}\approx 10^{-4}{M}_{\odot }{\mathrm{yr}}^{-1}$atR = (0.4–14) × 10¹⁵cm (forv_w = 25 km s⁻¹), our study points at a complex CSM density structure with asymmetries and clumps. The inferred densities are ≈10–100 times those of typical red supergiants, indicating an extreme mass-loss phase of the progenitor in the ≈200 yr preceding core collapse, which leads to the most X-ray luminous Type II SN and the one with the most delayed emergence of radio emission. These results add to the picture of the complex mass-loss history of massive stars on the verge of collapse and demonstrate the need for panchromatic campaigns to fully map their intricate environments. 

Abstract The end of supernova remnant (SNR) evolution is characterized by a so-called “radiative” stage, in which efficient cooling of the hot bubble inside the forward shock slows expansion, leading to eventual shock breakup. Understanding SNR evolution at this stage is vital for predicting feedback in galaxies, since SNRs are expected to deposit their energy and momentum into the interstellar medium at the ends of their lives. A key prediction of SNR evolutionary models is the formation at the onset of the radiative stage of a cold, dense shell behind the forward shock. However, searches for these shells via their neutral hydrogen emission have had limited success. We instead introduce an independent observational signal of shell formation arising from the interaction between nonthermal particles accelerated by the SNR forward shock (cosmic rays) and the dense shell. Using a semi-analytic model of particle acceleration based on state-of-the-art simulations coupled with a high-resolution hydrodynamic model of SNR evolution, we predict the nonthermal emission that arises from this interaction. We demonstrate that the onset of the radiative stage leads to nonthermal signatures from radio to gamma rays, including radio and gamma-ray brightening by nearly 2 orders of magnitude. Such a signature may be detectable with current instruments, and will be resolvable with the next generation of gamma-ray telescopes (namely, the Cherenkov Telescope Array). 

Abstract GeV and TeV emission from the forward shocks of supernova remnants (SNRs) indicates that they are capable particle accelerators, making them promising sources of Galactic cosmic rays (CRs). However, it remains uncertain whether thisγ-ray emission arises primarily from the decay of neutral pions produced by very-high-energy hadrons, or from inverse-Compton and/or bremsstrahlung emission from relativistic leptons. By applying a semi-analytic approach to non-linear diffusive shock acceleration, and calculating the particle and photon spectra produced in different environments, we parameterize the relative strength of hadronic and leptonic emission. We show that even if CR acceleration is likely to occur in all SNRs, the observed photon spectra may primarily reflect the environment surrounding the SNR: the emission is expected to look hadronic unless the ambient density is particularly low (with proton number density ≲0.1 cm⁻³) or the photon background is enhanced with respect to average Galactic values (with radiation energy densityu_rad≳ 10 eV cm⁻³). We introduce a hadronicity parameter to characterize how hadronic or leptonic we expect a source to look based on its environment, which can be used to guide the interpretation of currentγ-ray observations and the detection of high-energy neutrinos from SNRs. 

Abstract In 2021 August, the Fermi Large Area Telescope, H.E.S.S., and MAGIC detected GeV and TeVγ-ray emission from an outburst of recurrent nova RS Ophiuchi. This detection represents the first very high-energyγ-rays observed from a nova, and it opens a new window to study particle acceleration. Both H.E.S.S. and MAGIC described the observedγ-rays as arising from a single, external shock. In this paper, we perform detailed, multi-zone modeling of RS Ophiuchi’s 2021 outburst, including a self-consistent prescription for particle acceleration and magnetic field amplification. We demonstrate that, contrary to previous work, a single shock cannot simultaneously explain RS Ophiuchi’s GeV and TeV emission, in particular the spectral shape and distinct light-curve peaks. Instead, we put forward a model involving multiple shocks that reproduces the observedγ-ray spectrum and temporal evolution. The simultaneous appearance of multiple distinct velocity components in the nova optical spectrum over the first several days of the outburst supports the presence of distinct shocks, which may arise either from the strong latitudinal dependence of the density of the external circumbinary medium (e.g., in the binary equatorial plane versus the poles) or due to internal collisions within the white dwarf ejecta (which power theγ-ray emission in classical novae). 

Abstract Galactic cosmic rays (CRs) are accelerated at the forward shocks of supernova remnants (SNRs) via diffusive shock acceleration (DSA), an efficient acceleration mechanism that predicts power-law energy distributions of CRs. However, observations of nonthermal SNR emission imply CR energy distributions that are generally steeper than E −2 , the standard DSA prediction. Recent results from kinetic hybrid simulations suggest that such steep spectra may arise from the drift of magnetic structures with respect to the thermal plasma downstream of the shock. Using a semi-analytic model of nonlinear DSA, we investigate the implications that these results have on the phenomenology of a wide range of SNRs. By accounting for the motion of magnetic structures in the downstream, we produce CR energy distributions that are substantially steeper than E −2 and consistent with observations. Our formalism reproduces both modestly steep spectra of Galactic SNRs (∝ E −2.2 ) and the very steep spectra of young radio supernovae (∝ E −3 ).